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Upcoming gravitational wave (GW) detectors might detect a stochastic background of GWs po- 
tentially arising from many possible sources, including bubble collisions from a strongly first-order 
electroweak phase transition. We investigate whether it is possible to connect, via a semi-analytical 
approximation to the tunneling rate of scalar fields with quartic potentials, the GW signal through 
detonations with the parameters entering the potential that drives the electroweak phase transition. 
To this end, we consider a finite temperature effective potential similar in form to the Higgs po- 
tential in the Standard Model (SM). In the context of a semi-analytic approximation to the three 
dimensional Euclidean action, we derive a general approximate form for the tunneling temperature 
£SJ ■ and the relevant GW parameters. We explore the GW signal across the parameter space describing 

— ' the potential which drives the phase transition. We comment on the potential detectability of a GW 

signal with future experiments, and physical relevance of the associated potential parameters in the 



context of theories which have effective potentials similar in form to that of the SM. In particular we 
consider singlet, triplet, higher dimensional operators, and top-flavor extensions to the Higgs sector 
of the SM. We find that the addition of a temperature independent cubic term in the potential, 
CI arising from a gauge singlet for instance, can greatly enhance the GW power. The other parameters 

have milder, but potentially noticeable, effects. 

CD I. INTRODUCTION 

Gravitational wave (GW) interferometers currently taking data, and under development or construction, have as 
a possible target a stochastic background of GWs of cosmological origin (for reviews, see[l-4], for recent results, sec 
[5]). Since GWs propagate freely through the Universe after being produced, their detection provides a powerful 
diagnostics for the physics of the Early Universe. Several mechanisms that might generate such GWs have been 
discussed, including quantum fluctuations during or shortly after inflation [6-8], cosmic strings [!)], cosmological 
. magnetic fields [10-13], plasma turbulence [14-17] and bubble wall collisions during first-order phase transitions 
[18-23]. 

Specifically, the space interferometer LISA, expected to fly, or to be nearing completion, over the next decade will 
have a sensitivity peak for GWs with a frequency between 10 -4 and 1 Hz [4]. Quite fortunately, this range corresponds 
to the frequency range expected today, after redshifting, from GWs produced at a temperature T ~ 100 GcV ~ Eew> 
where the latter symbol indicates the electroweak scale. Intriguingly enough, the production of GWs at T ~ £"ew is 
indeed expected if the electroweak phase transition is strongly first order. In turn, this is a necessary condition for 
the success of scenarios where the baryon asymmetry of the Universe is produced at the electroweak phase transition 
(electroweak baryogenesis, EWB; for a review, see [24]). However, there is a generic tension between the requirement 
of a small bubble wall velocity needed in the context of electroweak baryogenesis models (see e.g. [25-27]) and the 
super-sonic bubble velocity values (detonation) under investigation here. Depending on the sources of net baryon 
number, super-sonic bubble wall velocities might be compatible with a baryon asymmetry produced at the EWPT, 
although this is generally not the case. We are at the dawn of the exciting period of exploring the electroweak scale 
with the possibility to connect such diverse experimental endeavors as GW detection and the Large Hadron Collider. 
These grand experimental enterprises may also help answer the question of how the baryon asymmetry arose in the 
Early Universe, which is all the more exciting. 

In its minimal setup, the strength of the electroweak phase transition in the Standard Model only depends only on 
the mass of the Higgs bosons (for a review, see [28]). For the Higgs mass range compatible with searches at LEP-II, 
non-perturbative lattice computations indicate that there is no phase transition at all, but rather a smooth crossover 
[291. No GW production is thus expected at the EW phase transition, if the electroweak sector corresponds to the 
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minimal Standard Model. In addition, if this is the case, EWB cannot be the mechanism underlying the production 
of the baryon asymmetry in the Universe. 

Fortunately, the strength of the EW phase transition is actually strongly model-dependent and parameter-dependent 
(for a review, see [30]), and numerous extensions of the Standard Model do predict a strongly first order EW phase 
transition. If this is further extended with the other necessary ingredients for EWB, the EW scale might indeed be 
responsible for generating the observed baryon-antibaryon asymmetry. 

If the EW phase transition is strongly first order, the universe finds itself trapped, at T < E^w, m a metastable EW 
unbroken phase where the vacuum state of the universe is a false vacuum (i.e. it is not the lowest energy state) as the 
universe cools down and its temperature decreases. A potential energy barrier exists between the false and the true 
(lowest energy) electroweak symmetry breaking vacuum. Quantum mechanical tunneling produces bubbles of true 
vacuum (broken EW phase), which then expand, collide, and combine to fill the universe with the true vacuum. GWs 
can be abundantly produced at the EW phase transition, primarily through bubble collisions [6-8, 17—20, 31]. This 
mechanism has been extensively investigated, analytically and numerically, and in relation to the effective potential 
that drives the transition itself. 

In particular, the possibility of a GW signal from the EW phase transition has attracted, not surprisingly, a great 
deal of attention. In view of the above mentioned large model-dependence, though, most of the existing literature 
has been devoted to either special cases and particular corners of parameter space (such as the light stop scenario in 
the context of the minimal supcrsymmctric extension of the Standard Model, MSSM; for a review, see [32]), and/or 
to accurate but numerical studies only [33, 34]. Alternatively, model independent results on the detectability of a 
GW signal from a strongly first order phase transition have been derived assuming a few relevant parameters for 
the GW dynamics could be computed from the effective potential that drives the particular phase transition under 
consideration (for recent related work, see e.g. [35-37]). 

The scope of the present study is to use semi-analytical results of the tunneling rate of scalar fields with quartic 
potentials [38] to predict analytically the strength of the EW phase transition and the GW signal in extensions of the 
Standard Model EW sector that can be characterized with the dynamics of a single order parameter, a scalar field 
4>. These models include simple generalizations of the SM effective potential, in appropriate dynamical regimes. The 
main result of this paper is a closed analytic formula for two parameters describing the GW signal as a function of 
parameters appearing in the effective potential of the scalar field driving the EW phase transition. 

The rest of this paper is organized as follows: in Section II we summarize the relevant physics and definitions of the 
quantities wc will be studying, in particular an approximation for the three dimensional Euclidean action, which is the 
basis for computing the finite temperature tunneling rate. Then, in Section III, we derive an approximation for the 
tunneling temperature and present exact and approximate formulas for the GW parameters for any effective potential 
similar in form to the SM case. Following that, in Section IV we constrain some parameters in order to reproduce the 
usual electroweak symmetry breaking pattern, and plot the effect of the various parameters in the potential on the 
parameters driving the GW signal. We examine the physical relevance of the parameter space and the detectability of 
the GW when varying various parameters beyond their SM values. In Section V we describe specific examples where 
our formalism can be applied in certain regimes, including SU (2) triplet and singlet extensions to the Higgs sector, 
top-flavor models, and models encompassing higher dimensional operators. 



In the context of a strongly first-order phase transition in the early universe, the basic problem of evaluating the 
GW signal amounts to calculating the tunneling rate (i.e. the decay probability) from the false vacuum state to the 
true vacuum state - in other words, the bubble nucleation rate per unit time per unit volume (for a review of tunneling 
in a finite temperature quantum field theory, see [28]). This is given in general by the expression: 



II. BACKGROUND 



r 



A{T)e~ s ^' T 



(HI) 



V 



where the factor A(T) ~ T 4 and Se3 is the three dimensional 1 Euclidean action 




(II.2) 



1 The symmetry for finite temperature bubble nucleation is 0(3), not 0(4) [39, 40]. 
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Typically Se3 is not calculated analytically, but semi- analytic expressions exist ([38]), which provide the basis for the 
rest of our analysis. We summarize these results below. 

As we point out in Sec. Ill, it is possible to derive an expression for the temperature at which tunneling is probable 
and the universe undergoes a phase transition. Now that regions of the universe can tunnel to the true vacuum, it is 
possible to produce large bubbles of true vacuum which expand, collide, and combine into larger bubbles. Expanding 
bubbles gain wall velocity and energy, but spherical symmetry does not allow any energy to be directly transferred 
into GWs. However, when more than two bubbles collide, this symmetry is broken and energy can be released into 
GWs 2 . Energy released into the universe can also be transferred to GWs through turbulence, but we neglect this 
sub-dominant contribution to the spectrum (for recent analyses, see [41-43]). 

Two quantities determine the GW spectrum when the phase transition proceeds through detonation (the bubble 
wall velocity is faster than the speed of sound in the plasma), a measures the energy density change by transitioning 
from the false to true vacuum, and /3 is the bubble nucleation rate per unit volume. The actual parametrization of 
the GW spectrum is summarized below. 

A. The Three Dimensional Euclidean Action 

The tunneling process is calculated from the three dimensional Euclidean action. Typically this is done numerically, 
however, there is a general, semi-analytic, approximate solution for quartic potentials of a single scalar field [38]. This 
is the key starting point for our calculations. 

Consider the potential for a scalar field cj) of the form 

V{4>) = A0 4 - a0 3 + b<t> 2 + af> + d (II.3) 

with A > to have the potential bounded from below, b > to have <f> = a minimum, and a > for the minimum 
to be for positive <f>. Without loss of generality, the false vacuum can be placed at the origin ((f> = 0, V = 0), and so 
d = c = 0. Here the above coefficients are typically not the same as in the zero temperature potential of the same 
theory, and may have temperature dependence. They are also all greater than zero. The three dimensional Euclidean 
action is approximately given by 

Sez = ^^(2 - Syi^^S + ^ + (3 3 S 3 ) (II .4) 

where S = 8Xb/a 2 , f3i = 8.2938, /?2 = —5.5330, and (3 3 = 0.8180. These parameters are the result of a numerical fit in 
the semi-analytic study 3 [38]. 

B. Gravitational Wave Parameters 

We are interested in the gravitational wave signal from electroweak symmetry breaking, which is mostly produced 
through bubble (as in which vacuum state) nucleation. The first part of calculating the GW parameters is to determine 
at what temperature bubble nucleation will be an energetically favored process. 

Assuming temperatures of O(100 GeV), the probability of a single bubble to be nucleated in a horizon volume to 
be ~ 0(1) is well approximated 4 in the Early Universe by 

S 3 (T t )/T t ~ 140 (II.5) 

where T t is the tunneling temperature (see e.g. [28]). This temperature will be between the critical temperature T c , 
where there is a degenerate minimum with <j> ^ (the only minimum at high T is at <fi = 0), and the destabilization 
temperature Td es t, where the minimum at <fi = is a local maximum: 

T dest <T t < T c . (II.6) 



2 Section IIB of [19] has a detailed discussion of the symmetry of two colliding bubbles. 

3 The absolute errors of the reduced action divided by the thin wall limit action are bounded to be less than 0.033 - see [38] for the full 
details. 

4 The exponential factor of the tunneling probability ensures that this approximation is valid for a broad range of temperature or energy 
scale. 
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An approximation is derived for T t in Section III. 

The vacuum energy (latent heat) density in this process is given by the standard statistical mechanics expression 



e t = -V(v(T),T)+T-^V(v(T),T) 

and the ratio between this and the radiation energy density is 

30e t 



(11.7) 



(II.; 



where g t is the number of relativistic degrees of freedom at T t . In a radiation dominated universe the parameter j3 is 
given by 



P d(S 3 (T)/T) 

-TT — J-t 



H t dT 



(II.9) 

T t 



a and (3/H t parametrize the GW spectrum [31], which is defined below, following [34]. 

We consider only a phase transition proceeding through detonation - the bubble wall velocity is faster than the 
speed of sound in the plasma. This also ensures that the thin wall approximation is valid, which was used in [34]. 
The wall velocity is [9] 



= v^ + V" + ^ (IL10) 

1 + a 

which increases with a, starting at the speed of sound in the plasma (1/V3) up to the speed of light. However, this 
may not always be a good assumption, depending on the exact theory for the electroweak phase transition; particle 
scattering with the bubble wall will affect Vb- For instance, [44] analyzed the MSSM stop contribution to the friction 
of the bubble wall in the plasma. This can greatly decrease Vb, down to about 0.05. On the other hand, one still 
expects some scaling with a (which can be quite large and we are not restricting ourselves to just the MSSM), which 
could counteract these effects. Since we do not make any assumptions on the underlying theory, it is difficult to 
say what value of u& will ultimately take, and we use the above equation for what follows. We are also restricting 
ourselves to the case of detonation, so we must have Vb greater than 1 / a/3- This is a non-trivial calculation from 
the effective potential and bubble dynamics, so for the purposes of this work we must assume that this holds. For a 
particular model, however, one must check that Vb is sufficiently large to apply the formulas we will present for the 
GW spectrum. 

How much of the vacuum energy is transferred to the bulk, rather than reheating the bubble, is given by an efficiency 
factor, again considering only the case of detonations and the same caveats above: 

*(«) = ^^77- f 0.715a + 4\/¥ I (H-H) 



1 + 0.715a 

For the bubble collision contribution the spectrum is parametrized (close to the peak frequency) as 

w = (IL12) 

°/Col* + a JCol* 

with the peak frequency fco\* arid peak amplitude f^coi* (as functions of the wall velocity). The exponents are in the 

range a £ [2.66, 2.82] and b £ [0.90, 1.19] with the case of large wall velocity, Vb ~ 1, having a sa 2.8 and b « 1.0 (see 
the numerical analysis of [34] for details). The spectrum observed today is found by redshifting: 

/c ol = 16.5 x 10- 3 mHz ( ^] ( ( T * r ) (JL) ^ (H.13) 

■ /Co1 \ f3 I \H t ) V 100 GeV / V 100/ v ' 



h 2 il Col = 1.67 x io- 5 o Col „ ( — 

\ 9t 



1/3 



H t \ 2 f a \ 2 /100 xl/3 
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with the functions /coi*//3 and A approximately [34] 

A = -MM, (11.15) 
0.42 + v 2 y ' 

^ = i.s-oZ + < (IL16) 

Although we use the results of a numerical analysis for the GW spectrum, there has also been recent analytic work 
on this subject (e.g. [36, 45]). 

To summarize, the three dimensional Euclidean action is used to calculate the tunneling rate to the true, symmetry 
breaking, minimum of the potential. An approximation for quartic potentials was given in [38], which we will employ 
in this study. The action is used to directly calculate the tunneling temperature, where this process can produce 
appreciable amounts of gravitational waves. Two parameters are calculated from the potential to give the spectrum 
observed today. 



III. T t ,a, AND p/H t FROM A GENERIC SM-LIKE POTENTIAL 

The potential we consider is a generic, quartic potential modeled after the SM Higgs's effective potential at high 
temperature. What follows is a quick overview of the relevant terms and quantities in the SM. The one loop, finite 
temperature, effective potential due to the gauge bosons and top quark, expanded at high temperature, is 

V(</>, T) = D(T 2 - T 2 SM )tf - ETf + (III.l) 



with the coefficients being 



D = 2ml ' ± m f + 2m ' « 0.17 (III.2) 
E = 2m w + m % w 9 6 x 10 -3 (nL3) 



n SM = ml ™* « (238.6 GeV) 2 (III.4) 
B = 73^4 ( 2 ™w +m%- 4m 4 ) « -4.6 x 10" 3 (III.5) 
X S m(T) = X SM ^ (a< log + m| log -gL - 4tt7 4 log J=L) (III.6) 

where logAs = logaf,— 3/2, log Af = loga/ — 3/2, a& = 167T 2 exp(3/2 — 2je), clf = tt 2 exp(3/2 — 2je) 7 and je ~ 0.5772 
is the Euler-Mascchcroni constant [46] (for details, see e.g. [28]). All the masses (Higgs, W, Z, and top) above refer 
to the usual zero temperature values, and v « 246 GeV is the Higgs vacuum expectation value (vev). The high 
temperature approximation comes from expanding the thermal bosonic and fcrmionic functions, Jb,f, appearing as 
the thermal contribution to the one-loop effective potential: 

|^Jb [m 2 (0)/T 2 ] , ^|pJ F [m 2 (</>)/T 2 ] , (III.7) 

where m is the field-dependent mass of the boson or fcrmion and we are working in units where k = 1 = h = c (so 
temperature is measured in energy units). The thermal functions are given by 



Jb,f [m 2 f3 2 ] = / dxa; 2 log [lTe-^ 2+m2/32 , (III.8) 
and have high temperature expansions 

4 2 9 / 9 \ 3/2 1 4 9 

MmVT 2 ) + -- w]oe — + ... (III.9) 

. o , 9, 77T 4 7T 2 771 2 1 TO 4 777 2 . 

MmVT*) « - + -- - -- A log _ + • • • (111.10) 
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Since T t > Tqsm, the destabilization temperature, the temperature scale for the phase transition is much larger than 
any of the particle masses; a high temperature expansion will be a good approximation in this regime. 
In terms of the definitions given in Section II we have 



A(T) = 



a{T) = ET 
b{T) = D(T 2 



Tqsm) 



5(T)- 8Xb - 2X sm(T)D(T 2 -T ( 

n 2 



2 

OSM 



E 2 T 2 



2D 

£2 



1 - 



Tqsm \ ' 



T 



Asm(T). 



We now consider an effective potential of the general form of the SM case 



v, 



eff\ 



- (ET - e) 



D(T 2 - T 2 ' 



(111.11) 
(111.12) 
(111.13) 

(in. 14) 



(111.15) 



where there is a new parameter, e, motivated by gauge singlets (see Sec. V A), e = in the SM. The potential minima 
(first derivative with respect to <p is 0) located at 



= 



3(ET - e) ± y/9(e - ET) 2 - 8D\(T)(T 2 - T 2 ) 
2A(T) ' 



(111.16) 



The parameters in the Euclidean action are (using a tilde here to denote the A of the formulae in Sec. II A, otherwise 
we moan the A of our general effective potential above) : 



ACT) 
4 

ET - 



5(T) = 



A(T) = 

o(T) = 
b(T) = D(T 2 - T 2 ) 

8\b _ 2D(T 2 - T 2 )A(T) 
~^ ~ {ET - e) 2 ' 



(111.17) 

(111.18) 
(111.19) 

(111.20) 



The critical temperature is (the sign of the square root is chosen so that the temperature is positive in the SM case): 

eE 



T- 



y/gAggg + (DX(T C ) - E 2 )fJ) 
E 2 - DX(T C ) 



(111.21) 



From here on the temperature dependence of A will be dropped, taking it as a free parameter in the theory. Alterna- 
tively, A(T) rj A(T ) or A(T) sa A(T C ), since A is slowly varying with temperature (logarithmic corrections). 

For calculating the tunneling temperature, we note the general behavior of Ses/T: it decreases rapidly as the 
temperature is lowered, from a singularity at T = T c . Physically, this is because the tunneling rate starts at zero when 
the minima are degenerate, and increases rapidly as the new minima becomes the global one. In the approximation, 
this can be seen by looking at the term (2 — S)~ 2 . At T c , 8(T C ) = 2, for this general effective potential. To 
approximate a solution for T t we expand near T c : T — > T c (\ — e). e will be very small for T t , as the sharp peak ensures 
that Ses/T = 140 very close to T c . Expanding 5 to lowest non- vanishing order in e, about e = 0, 



4D\(ET 2 - eT c )T c 
+ (e — ET C ) 3 6 



Ft, 



(111.22) 



where the first term is due to substituting in (III. 21). As e — J- 0, we recover that 8 — > 2. Now the relevant quantity for 
the tunneling temperature, Ses/T, can be approximated by again having T — > T c (l — e) and expanding about e = 0. 
Defining the prefactors of Se3 as 



G = 



64V^tt 



the resulting lowest order expression is 

Ses/T: 



81A 3 / 2 ' 

2G(P 1 + 2I3 2 + A/3 3 )(ET C 
F 2 T C 



e) 1 



(111.23) 



(111.24) 



7 



As expected, there is a singularity as e — !• 0, with the same power as the divergent piece in the original expression. 
Finally, we arrive at an estimate for the tunneling temperature: 



( * V 140F^ (IIL25) 

T t « T c (l - e). (111.26) 

A rough error estimate comes from comparing this value of T t to the value obtained numerically using the first 
approximation, (II. 4). Throughout most of the parameter space that we analyze in the next section, the difference 
between the two values averages to be less than 0.1% (most of the space is much lower even). Some regions of — e, 
however, can have an average difference of 30%. 

The true minimum of the potential is <\> with the positive square root in (III. 16) which gives an expression for the 
minimum of the potential as a function of temperature. This yields an exact expression for a (from the formulae in 
Sec. II B): 

a = 15( i6~2g t A^° < 9e3 + e2{9ETt ~ 60 + ^ T t^ET t + - 4Z?A [^ET t {T 2 - 3T 2 ) 

+ (T 2 + T 2 )(\ + e [-45£ 2 T 2 + 12L>A(T 2 + T 2 ) - 3ET t £ + £ 2 ] } , (111.27) 



where £ = y/9{e - ET t ) 2 - 8L>A(T t 2 - T 2 ). @/H t can be computed directly from (II.4) and the definitions in Sec. Ill: 



PIH t = 1024ttJ 2 V2T ^_ r_ 32DXT ,( ET 2 _ eT )( ET _ e )h 2 j(8, 64) 

' 81T t (e-ET t ) w VDXhk^ X V ° A ' V ' 
-ET t h 2 kj{8, 64) + (ET t - e)h 2 kj{8, 64) + T t (eT t - ET$)hkj(8, 64) 

+4T t (ET$ - eT t )(e - ET t ) [(e - ET t f - ADXh] hj(W, 192)} , (111.28) 

where 

h = T 2 - T 2 (111.29) 

j(x, y) = [ET t - e) 4 A + xD\(ET t - efhfo + yD 2 \ 2 h 2 p 3 (111.30) 

A much simpler expression comes from using the approximation in (III. 24) (noting that = — yrg^ near T t ): 

a /EJ ^, 4G(/3! + 2/3 2 + 4/3 3 )(igr c - e) (1 - e) 

P/-«t ~ ^7 — ■ (111.32) 

Again, a rough error estimate comes from comparing this approximation versus computing f3/H t directly from (II. 4). 
Over the parameter space we consider, the difference typically averages to be less than 10%. Again, part of the — e 
region has much larger errors, but, on average, well within an order of magnitude. These errors in T t and P / H t might 
potentially affect the peak frequency and overall power, to an extent that quantitatively depends on the size of these 
errors. It is difficult to estimate analytically the effect on the GW spectrum, but if the errors are not very large, they 
are possibly only important in borderline detection cases. 

Using the above approximation for the tunneling temperature provides an approximation for the GW parameters 
for any theory with a SM-like effective potential. The GW spectrum is then fully specified, using the equations at the 
end of Sec. II B. 



A. Parameter Constraints 



We need to enforce that the potential of (III. 15) (where we will work at T = here) correctly describes electroweak 
symmetry breaking. Note that the potential is typically not a tree level potential, even at T = (in the SM there are 
still the one-loop effects), and using the SM parameters will give differing results (by a few percent) from tree level 
calculations. The first condition is that <p a t the electroweak breaking minimum, (III. 16) at T = with the positive 
root, is the usual vev of the Higgs field: 



< >= - 3e+ V /9 e ; + 8mT ° = „ « 246 GeV, (111.33) 
2A 
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where A is also at T — 0. This point must be a stable minimum, hence 



d 2 V eff ^,T = 0) 



-!- ( 9e 2 + 8L>AT 2 - 3 ev /9e 2 + 8DAT 2 ) > 0. (111.34) 



Finally, we will also restrict the parameters to have a Higgs mass above current limits 

2 _ d 2 V eff (<f> 7 T = 0) 



> (114GeV) 2 (111.35) 

(p—V 

The first constraint, (111.33), can be solved to give an expression for Tq: 

n = v -^^ (in.36) 

which is rather similar to the SM form. Requiring that To and v are greater than 0, along with the original constraints 
on the coefficients of the potential (A, D, (ET—e) > 0, see Section II A and [38]), satisfies the second constraint, (III. 34), 
if e > or if e < and A ^ — Finally, we can use the the last equation for m^, (III. 35), to solve for A. In order to 
satisfy all the constraints we are left with the following solutions: 



m h —3ev 



5 4^( m l ~ 9ew _ V m h ~ l&evml + 9e 2 v 2 ), 
if ( !! ^' =T { 3v + V<K + ^ 2 )) < e < 0, (111.37) 
m Jf^. : if <e<^. 

2v z 1 — 3v 

For the SM case of e = we have the usual relation of m 2 = 2 Aw 2 . In order for the theory to be perturbative we need 
A < 1, so for the SM we use 

0.11 < A < 1, (111.38) 

where the lower bound is from the Higgs mass limit. This roughly corresponds to 

115 GcV < m h < 348 GcV. (111.39) 

When e ^ the value of e and nih in the range given above set the value of A. When e < both solutions for A arc 
used, which gives the greatest range of values for both A and e. 

Other constraints follow from requiring T c to be real and positive, which can be deduced simply from (III. 21). For 
instance, for T c to be real, 

e 2 > (E 2 — DX)Tq, (111.40) 

while having T c > can reduce to a simple constraint depending on which parameters are being varied together, as 
well as the sign of the square-root. 

It is also important to note divergent features of the above expressions. First, in (III. 25), e will be complex if 
e > ET C > 0, since the other terms are all positive (D and E are positive in general, as they depend only on 
masses-squared). In this case T t and the GW parameters will also be complex, so we limit 

e < ET C . (111.41) 

Since this arises in our approximation, it is possible that higher order terms will resolve this behavior. However, if 
e > ET C , then the sign of the cubic term in the potential changes (becoming positive). There would no longer be any 
potential energy barrier, and the phase transition could not be first-order. Therefore, it is consistent for this analysis 
to limit e by the inequality above. 

In several of the equations in Section III the term E 2 — DX appears. Occurring in the numerator of eq. (III. 21) 
this appears in the reality constraint above. Appearing in the denominator, however, when it equals zero it is a pole 
in the (exact) expression for T c . This causes e to diverge, and thus T t will pass through zero. This carries over to a 
pole in a. Clearly, as e approaches one, our approximation will break down. However, fine-tuning to be near this pole 
could provide models with a very strong phase transition (large a). Below, we only consider values of E 2 up to 90% 
of DX. 
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FIG. 1: In each plot the labels for the x-axis denotes the varied parameter in blue and red (dashed), respectively. The green 
line is the (constant) SM value. The ranges for the varied parameters were chosen arbitrarily, to show the overall behavior, e 
is in GeV, while the other parameters are dimensionless. Note the pole as E is varied, as described. On the left-hand side, the 
blue line (varying e) closely follows the green one (SM value), for the range shown. 



IV. THE PARAMETER SPACE 



Now that we have general expressions for the GW parameters in terms of the coefficients of our potential, we can 
look at the parameter space. Throughout this analysis the parameters we are not varying are set to their SM values 
(with \sm = 0.3, and \sm(T)\ t _ t when not varied). Also, we set g t — 100 throughout. This can vary greatly 
between models, but here we are looking at the parameter space in general, without assumptions on what theory is 
being used. 

In the plots of the a — ft / H t plane both A and one other parameter are varied. A ranges as described above. The 
other varied parameter ranges from an order of magnitude less than its SM value to an order of magnitude larger 
for the darker blue region, and a further order of magnitude larger and smaller for the lighter blue region, taking 
into account the above constraints. Note that for negative e, we also need to make sure that e is real and less than 
one, as this is the only parameter (in the ranges we considered) that could cause these problems for e. The red line 
represents the SM with only A varied, e is less than ET C or m^/3i>, whichever is smallest. LISA and BBO sensitivity 
regions were computed using rough estimates of the sensitivity curves, as in [34, 35, 43, 47], and finding the region 
in the a-/3/H plane where the computed GW spectrum is above the sensitivity. The LISA spectrum used is the 
approximate instrument strain sensitivity, while the BBO spectrum is the correlated extension (BBO Corr), which 
uses data analysis of two LISA-like detectors to improve sensitivity. LIGO is not sensitive in this region to GWs. 
The sensitivity regions shown flatten (do not reach into higher 0) for large a, as the LISA and BBO sensitivity drops 
considerably outside of their most sensitive frequency. Thus, even with increasing a, larger j3 means a greater peak 
frequency, pushing it out of the experimental range. 

We note that very small Vb, due to viscosity of the plasma, could potentially affect detection with LISA or BBO, 
due to a change in the peak energy density and the peak frequency of the GWs. A small Vb = 0.05, as in [44], does 
not change the overall overlap of the experimental sensitivities with the parameter space in the figures (using this 
Vb in calculating the experimental sensitivity). However, this could be a large effect for borderline cases. Without 
specifying the specifics of a theory, it is difficult to explore the issue of calculating Vb precisely. Furthermore, our 
analysis considers only detonations, where is at least ten times greater. Small requires considering the case of 
deflagrations in detail, and is outside the scope of this work (for recent work, see [36, 48]). This is another direction 
to be explored not only in the context of GW detection, but also for baryogenesis, as Vb affects the length of the phase 
transition. 
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FIG. 2: A plot of the a — {3/ H plane with A and D varying. The red line is the SM with A in the allowed range. The darker blue 
region is an order of magnitude greater and less than the SM value for D. The lighter blue region is a further order of magnitude 
greater and less. BBO and LISA rough sensitivity regions are shown as the lighter and darker red shading, respectively. 




10 10 10 7 10 4 0.1 100 10 5 



FIG. 3: A plot of the a — j3/H plane with A and E varying. The red line is the SM with A in the allowed range. The darker blue 
region is an order of magnitude greater and less than the SM value for E. The lighter blue region is a further order of magnitude 
greater and less. BBO and LISA rough sensitivity regions are shown as the lighter and darker red shading, respectively. 
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FIG. 4: A plot of the a — /3/H plane with A and e varying. The dark blue region is for e > 0, while the light blue shading 
indicates the first solution for A when e < 0, and the green the second solution. BBO and LISA rough sensitivity regions are 
shown as the lighter and darker red shading, respectively. The tail of the green region, extending to very large a, may have 
errors in the calculation of P/Ht due to the approximate formula. However, the points remain in the LISA region. 



Although varying D or E docs not greatly enhance a, their effect is noticeable. Both also tend to have large values 
of /3, but cover a very large range. However, D, at least for this range of variation, is not close to the sensitivity 
regions of LISA or BBO. E has regions which are much closer to LISA or BBO, but we do not predict any overlap, 
e < 0, on the other hand, can greatly increase a, as well as have a much smaller (3. Here there is the greatest 
parameter space which may be observable. There is considerable overlap between the various potential parameters 
and the corresponding values of the GW parameters, and thus a particular GW spectrum does not generally point to 
a value of just one of these potential parameters. 

V. MODELS 

In this section we investigate a few models which naturally fit into the effective potential analyzed above. Of 
primary interest are models with gauge singlets, which we argue provide a very good motivation for the parameter e. 
Since this parameter was shown to have the largest effect in producing a strongly first-order phase transition, these 
models could predict an observable GW signal in the near future. There are many models with singlets, and in general 
it is not possible to derive an analytic expression for the effective potential. We will look at the simplest such model, 
the SM with the addition of one (real) gauge singlet field. 

A. SM + Gauge Singlet 

While the motivation for the overall form of the effective potential we consider is due to the SM, the additional 
term, e, is motivated by gauge singlets. Gauge singlets are present in many models, such as the Next-to-MSSM 
(NMSSM). Singlets in the context of the electroweak phase transition and baryogenesis have been studied for some 
time (e.g. [49-51], and more recently, e.g. [52-54]), including the effect of trilincar couplings driving a strongly first- 
order phase transition without even considering the thermal term (see, for instance, [55, 56]). Recently, [54] has also 
explored the connection of singlets with possible explanations for recent astrophysical signals. In this section we will 
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give a brief decoupling (in the sense of taking the singlet to be much heavier than the Higgs) analysis, motivating 
the inclusion of e in the effective potential, and showing that these classes of models naturally fit into our framework 
above. 

The procedure we will follow is to use the tree-level, zero-temperature mass eigenstates as a basis, taking the mostly 
singlet state to be much heavier. In addition, we will want the mixing of the states to be very small. The goal is 
to show that there is a tree-level, cubic, self-coupling term of the lighter (mostly Higgs) state. We will consider two 
different limiting cases (in terms of an expansion parameter, defined below) and show that we can always have a 
very small mixing angle, a very heavy singlet state, and a light Higgs state with sufficiently large cubic coupling. If 
we consider SM values for all other parameters of the effective potential, than we need e < — 20 GcV (see previous 
section) for a strongly first-order phase transition. 

Following the notation of [53] , the tree level scalar potential of the SM Higgs sector with an additional gauge singlet 
field is the sum of 



VsM 



a i 



2 (fftm 



A (H*H) 
a 2 



V HS = f(^H)S+f(H^HyS 2 
Vs = ^S 2 + |s 3 + h -±S* 



(V.l) 
(V.2) 

(V.3) 



with H the usual SM SU (2) l scalar doublet and S the gauge singlet (real) scalar field. The vevs are defined to be 
Vq/\/2 for H and xo for S. The fields for the fluctuations about the vevs are h and s, defined as H = (vq + K)/y/2 and 
S = xq + s. For simplicity, we'll now take H to be real. The minimization conditions (for both fields, with xo =/= 0) 
can be used to eliminate the mass parameters: 
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(V.4) 
(V.5) 



While for xq = the conditions enforce a\ = and /i 2 
The mass matrix has the following elements: 



XsVq, as in the SM. 
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(V.6) 
(V.7) 
(V.8) 



The mass eigenstates are defined as 



hi = (sin#)s + (cos 6)h 
hi = (cos8)s — (sin6)h 



(V.9) 
(V.10) 



with the mixing angle 6 as 



tan# = 



1 + + 



where y 



2 



(V.ll) 



Here, |cos^| > l/\/2, and so h\ is the state with the largest 5f/(2)-like component (and /12 has the largest singlet 
component). The terms singlet state, hi, and heavier state will be used interchangeably (and likewise for the Higgs 
state). Inverting the above states, we have the original fields, expanded about the minimum, in terms of these mass 
eigenstates: 



H = vq + (cos#)/h — (sin 8)h 2 
S = xq + (sm6)hi + (cos9)h2- 



(V.12) 
(V.13) 



The mass eigenvalues are 



m 



2 

1,2 



A + n 



± 



(V.14) 
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with the upper (lower) sign for mi (771,2). 

We consider the decoupling limit where the singlet is very heavy (i.e. the state h 2 ) and study the cubic term of the 
effective hi potential. The lighter state, hi, should roughly be in the mass range allowed for the SM Higgs. Although 
we allow mixing between the states (so the singlet will not be completely physically decoupled), we want to keep the 
mixing angle small; interactions between hi and h 2 will be highly suppressed after integrating h 2 out, by both the 
small mixing angle and large mass scale of h 2 . However, since higher dimensional operators are generated, this could 
contribute to enhancing E (and thus a), as discussed in [57]. We will also consider xq an essentially free parameter, 
determined by physics at the higher, singlet scale. We will consider xq both much larger or smaller than vq, and use 
the ratio as an expansion parameter. 

Since we consider xq set by other dynamics, S will be expanded as in eq. (V.12), while we will drop the Vo in 
expanding H b . The reason for this is that we want to consider these results as leading terms in the finite-temperature 
effective potential; the proper degree of freedom before the phase transition is the unshiftcd field (there is no vev yet). 
Clearly, the mass eigenstates above are for the tree-level, zero-temperature, shifted fields. However, we will consider 
this as simply a change of basis; the coefficients of the expansion are just the proper values for the mass eigenstate 
basis at tree level (and zero temperature). 

Considering the limit of xq 3> vq, we expand to lowest order in the small parameter u = vq/xq. We have that 



cos# w 1 and sin 6* w —^§-u. Also, to lowest order in u, 



2b 4 
Mo 



ir 



(V.16) 



Note that as we take u very small, the heavy mass becomes arbitrarily large. Writing the fields in the potential in 
terms of the mass eigenstates and working to lowest order in u, the coefficient of h\ is given by 

a\v Q a 2 ( a 2 b 3 \ , V17 , 
— ai I u. ( V .17) 



2b 4 2b A \ 26 4 

Some fine tuning will be needed to have appropriate masses and cubic coupling. The values, for example, of a 2 = 
0.35, 64 = 0.3, A = 0.4, has a light mass of about 190 GeV and a cubic coupling of about —50 GeV (in the limit of 
u — > 0). The massive parameter ai could also be tuned as 1/u, for instance. In this case the series expansions also 
change and there are large regions of parameter space where all the constraints can be met. 
In the small Xq limit, to lowest order in w = Xq/vq, 

6k2w (V.18) 

m\ w AXvl + aiv w (V.19) 
2 ai^o 



4w 



(V.20) 



Again, the heavy state mass is arbitrarily large as we take w very small. The cubic coefficient of hi, in this limit, is 
simply aiw. In this case, we can have m 2 large and a sufficient cubic coupling by having ai negative and (possibly 
unnaturally) large. However, the second term in the expression for m\ can constrain ai (or require tuning with A). 
For instance, the light state has mass about 177 GeV when A = 0.3 and ai = —20/w (the cubic coupling is —20 GeV), 
with w arbitrary. Note that these expansions (to lowest order) do not change for ai = a/w. 

In both limiting cases then, a cubic coupling in the effective potential of the light state appears. This coupling can 
be of the right size to drive a strongly first-order phase transition, while still having a small mixing angle between the 
states, and appropriate heavy/light masses for the fields. It is also possible (for instance, [53]) to approximate the 
two-field effective potential in the same form as we analyzed (with the field being the Higgs, not the singlet). Viewed 
this way, the coefficients are more general functions of the singlet field. In general, however, this potential still must 
be numerically analyzed, due to the dynamics of the Higgs-singlet interactions. The decoupling analysis shown above 
is simply one limit of a general analysis. It is meant to motivate e and show that, even in this extreme limit, a singlet 
can have a large effect on the phase transition 6 . 



5 If we include vo, the only change is the addition of the expected cubic coupling of the shifted Higgs field in the SM. 

6 The singlet may also have other effects on the overall potential, including from its own finite temperature potential. A full analysis, of 
course, must account for effects besides generating e. 
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B. SM + SU(2) L Triplet 

An example of a model where the quartic coupling in the tree-level effective potential is suppressed for a fixed value 
of the SM-like Higgs is given by models with an additional SU(2)l triplet in the scalar sector. First considered in [58], 
this extension to the SM Higgs sector has several significant phcnomcnological implications, including a dark matter 
candidate [59, 60] and providing a natural framework for a Type II seesaw mechanism for generating non-vanishing 
neutrino masses [61]. We postpone an exhaustive study of this scenario to a future paper [62], but we point out here 
this scenario as one where the results of the present analysis apply. 

Indicating with S the neutral component of the SU (2) l triplet A, the tree level neutral CP-even part of the Higgs 
potential, neglecting interaction terms with two or four A's, reads: 

V(H, 6) = -m 2 H H 2 + ^H 4 + M 2 A S 2 - 2 f ,6H 2 , (V.21) 

where Ma is the mass term associated to A and the last term stems from the term, in the scalar potential, 

HH T ia 2 £JH + h.c. (V.22) 

(see e.g. [61]). Imposing the minimization condition to the tree level potential and singling out the field-space trajectory 
along which the minimum of the (tree-level) potential is found allows us to express 5 as a function of H: 



— = 5=J^v 
85 M 2 



H- 5 =T7tH 2 - (V.23) 



Along the locus of minima in the 8 direction, the effective potential is now a function of a single field, H: 

V min (H) = -m 2 H H 2 + - |1) H\ (V.24) 

This is essentially the same as in the SM case, but with the parameter change 

Asm Asm 



_ A T . (V.25) 



This model can therefore be cast into the form of the effective potential of the SM we consider in this study, with 
the above change to the effective quartic coupling A. The extra triplet scalar therefore automatically enhances the 
strength of the phase transition, since, for a given SM Higgs mass and at a given tri-scalar coupling E, a smaller A 
correspond to larger values of a (which we can think of as the "strength" of the phase transition). 



C. Other Models 



Here wc will only briefly mention a few other models which could be analyzed in our framework. First is the 
analysis of [63]. The authors study baryogencsis using the MSSM to generate a low energy effective theory, and in 
particular also study the strength of the electroweak phase transition. Given certain constraints on the parameters, 
the potential minimization reduces to the one-dimensional case, and the phase transition strength (i.e. the value of 
a) can be greatly enhanced. It is enhanced through an increase of the parameter E, which can be about an order of 
magnitude larger than in the SM. Given SM values for the other parameters in the effective potential, this alone, for 
a light Higgs, drives a to be larger than one. The new expressions for E can also be used to investigate more closely 
the parameter space that produces large a. 

In a similar vein, [57] considers dimension-six Higgs operators, arising, for instance, from integrating out a heavy 
singlet. This can also produce a first-order phase transition. Another effect of this operator is to alter the Higgs self- 
couplings from the SM. The cubic and quartic self-couplings are altered, and this in turn would alter the couplings 
appearing in the effective potential. Again, this is easily incorporated in our analysis. 

The "topflavor" model also has a phase transition which fits into our framework. In this model there are separate 
SU(2ys for the third versus other generations. In [64], the earlier phase transition, from SU(2)i x SU{2)2 to the 
SM SU(2)l x U(1)y, was analyzed in the context of baryogenesis. This phase transition has a scalar field as the 
order parameter, which has a quartic (tree-level) potential. The one-loop, finitc-tcmpcrature, effective potential can 
be derived, and it has the same functional structure as the SM. The high temperature expansion will then be the 
same form as the general effective potential we analyzed. By using the constants of the topflavor model, which has a 
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strongly first order phase transition, the GW parameters can be found through the above results. Note that in this 
case some of the constraints on the parameters are not applicable, since they refer to specific electroweak constraints. 
Baryogcnesis in this model requires that the electroweak phase transition is not strongly first order, and therefore any 
GW produced comes from this earlier phase transition. 

Finally, we note that the form of the potential we used is really quite general. The potential mimics the SM form, 
which we might reasonably expect to be a good approximation to any low energy, effective theory. Also, besides 
e, the form of the potential and the constants are all derived from the high-temperature expansion of the general 
one-loop, finite-temperature, effective potential of a theory with gauge bosons and fermions coupled to a scalar field. 
As we have shown above, the addition of gauge singlets motivates the inclusion of the additional parameter. Thus we 
expect this potential to arise very generically for any field theory model. The obvious exception is for multiple fields 
controlling the phase transition (so that all but one cannot simply be integrated out). In that case, it may still be 
possible to use this potential, by approximating field configurations or by taking certain limits. 

VI. CONCLUSIONS 

Upcoming experiments may soon observe the first GW signals, including those from early universe processes. Such 
experiments will deliver the very exciting prospect of observing GWs from the electroweak phase transition, if it is 
strongly first-order. While there have been many studies on this topic, they largely involve numerical computations 
of the phase transition temperature and GW parameters. In this work we have presented an analysis of a generic 
effective potential which is motivated by the SM Higgs potential and an additional contribution from a gauge singlet 
(or other possible origins) in the form of a temperature-independent cubic coupling. 

By approximating the tunneling temperature to be very close to the critical temperature (based on the general form 
of the action) , we have derived an expression for the tunneling temperature based on the parameters of the effective 
potential. Using this result, the GW parameters also have expressions depending on the parameters of the potential. 

Once these expressions were found, the parameter space of the potential was explored starting from the SM values. 
While all the parameters can have noticeable effects, e, motivated from singlet models, easily has the greatest effect. 
Following this, we showed how this parameter arises from a decoupling analysis of a simple SM plus singlet model. 

Other models can also be analyzed with our effective potential. The addition of an additional SU (2) triplet affects 
A, while higher dimensional operators can also affect E. Increasing E by an order of magnitude from its SM value 
could also produce a strongly first-order phase transition. The topflavor model also has an effective potential of the 
same form as the SM for its earlier phase transition, which can also be analyzed through this work. In this case, the 
GW signal would not be from the (later) electroweak phase transition. 

Finally, although the effective potential studied has clear origins in the SM and singlet extensions, the potential 
is indeed very general. A low energy, effective theory attempting to model the electroweak phase transition will 
likely closely model this form of the potential. Additionally, the potential is a high-temperature approximation to a 
one-loop, finite temperature, effective potential of a scalar field with gauge bosons, fermions, and singlets. So in this 
sense its form is also generic and expected to be common in a phase transition with a scalar field order parameter. 

The potential applications of the results presented include applying it to other models for the electroweak (or other) 
phase transition. By approximating or otherwise finding a way to put an effective potential in the form we analyzed, 
the tunneling temperature and GW spectrum is now easily obtained through the above formulae. Besides being 
significant itself, the detection of a stochastic background of GWs may provide insight into the specifics of models of 
the electroweak phase transition, and the fundamental mechanism underlying the generation of the baryon asymmetry 
of the universe. 
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